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Abstract 


Fluorescence  emitted  by  an  atom  near  a  metal  surface  and  the  coupling  of  the 
surface  plasmon  field  of  a  dielectric  substrate  to  a  molecular  electronic 
transition  are  studied.  Explicit  expressions  for  the  atomic  and  molecular 
lifetimes  are  derived.  It  is  shown  how  the  (classical)  phase  shift  of  a  plane 
wave  upon  reflection  at  the  surface  is  responsible  for  the  alteration  of  atomic 
lifetimes.  Subsequently,  it  is  demonstrated  that  the  dipole  direction  of  an  atom 
can  be  fixed  by  illumination  of  the  system  with  a  polarized  light  source.  For 
molecular  transitions,  the  surface-plasmon  and  the  surface -roughness  contribution 
to  the  decay  constants  are  included.  In  a  comparison  between  the  image  theory 
and  the  energy-transfer  theory,  it  appears  that  the  former  can  be  rather 


I.  INTRODUCTION 


An  atom  in  empty  space  decays  spontaneously  (accompanied  by  the  emission  of 
a  fluorescent  photon),  with  a  certain  lifetime,  which  is  brought  about  by  the 
coupling  of  the  atomic  (transition)  dipole  moment  to  the  travelling  plane  wave 
modes  of  the  electromagnetic  vacuum  field.  Early  attempts  to  comprehend  the 
process  of  spontaneous  emission  concentrated  on  the  prediction  of  the  Einstein  Aq 
coefficient,  which  equals  the  inverse  lifetime  of  an  excited  state,  and  it  was 
one  of  the  first  victories  of  quantum  mechanics  by  Dirac  in  1927  to  establish 
unambiguously  its  value  as 
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(1.1) 


Here,  <j  ||y||j  >  is  the  reduced  matrix  element  of  the  atomic  dipole  operator  y 
eg  ~ 

between  an  excited  level  |j  m  >  and  the  ground  level  |j  m  >,  which  have  a  level 

1  e  e  8  8 

separation  Kw.  Result  (1.1)  follows  from  the  Golden  Rule  and  is  an  essentially 

quantum  mechanical  expression  (because  of  the  appearance  of  ft). 

1  2 

Earlier  this  year  *  it  was  demonstrated  for  the  first  time  that  it  is 
practically  feasible  to  alter  atomic  lifetimes  for  optical  transitions  through  a 
(semi)  confinement  of  the  atom.  This  was  brought  about  by  passing  an  atomic  beam 
through  a  set  of  parallel  mirrors  or  a  confocal  resonator.  In  this  paper  we 
report  recent  developments  in  the  theoretical  under standing  of  atomic 
fluorescence  near  a  metal  surface. 

Another  technique  to  change  the  lifetime  of  an  electronic  transition 
pertains  to  the  deposition  of  a  molecular  monolayer  on  a  dielectric  bulk 
material.  Then  the  separation  between  the  optically-active  compound  and  the 
passive  substrate  can  be  ^extremely  small,  which  implies  that  the  lifetime 


modification  is  brought  about  mainly  by  coupling  to  tha  short -rang*  avanascant 
plasmon  fiald.  Sinca  tha  decaying  system  is  so  close  to  tha  surface,  it  becomes 
inevitable  to  take  the  details  of  the  surface  structure  (roughness)  into 
consideration. 


II.  ATOMIC  DECAY 

The  plane  z  ■  0  separates  the  vacuum  z  >  0  from  a  substrate  in  s  <  0,  which 
is  assumed  to  be  fairly  represented  by  its  dielectric  constant  c(<i>).  An  atom  is 
positioned  on  the  z-azis  at  de  with  d  >  0.  Since  every  radiation  field  in  the 
region  above  the  surface  can  be  expanded  in  plane  waves,  it  should  be  possible  to 
express  any  surface  effect  in  an  observable  quantity  in  terms  of  the  plane-wave 
Fresnel  reflection  coefficients.  These  are  for  s(urface)  and  p(lane)  polarized 
waves 
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where  u  signifies  the  cosine  of  the  angle  of  incidence.  Then  the  mode  expansion 
of  the  electric  field  in  the  region  z  >  0,  becomes 


E(r,t) 
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with  F.  (r)  the  sum  of  an  incident  and  a  reflected  wave,  and  normalized  in  such  a 
~zo  - 
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way  that  avary  photon  carrias  an  anargy  Tha  aquation  of  motion  for  tha 

danaity  oparator  of  tha  aystam  raada 

i*^  -  tHa  ♦  Hr  «•  Rt  ,  pj  .  (2.4) 

whara  a,  r  and  i  indicata  atomic,  radiation  and  intaraction,  reapactivaly. 

Transparent  intarpratations  can  ba  obtainad  if  wa  raatrlct  our  attantion  to 
(idaalisad)  metals,  defined  by  c(«i)<0.  From  Eqs.  (2.1)  and  (2.2)  it  then  follows 
that 

|Ro(w,u)|  -  1  ,  (2.5) 

and  hence  tha  metal  is  perfectly  reflecting.  Than  wa  can  write 

214  (u,u) 

R  (u,u)  •  a  ,  (2.6) 

o 

and  tha  raal-valuad  phase  shifts  $g  for  a  and  p  waves  follow  directly  from  a 
comparison  with  Eqs.  (2.1)  and  (2.2).  For  a  perfectly-conducting  metal  wa  have 

♦.  *  *• 

In  writing  for  a  magnetic  substata  of  tha  upper  level,  wa  implicitly 

refer  to  a  direction  in  configuration  space,  e.g.,  tha  quantisation  axis.  In  our 
geometry  it  is  most  advantageous  to  taka  this  direction  to  coincide  with  tha 
normal  to  tha  surface.  With  this  convention  the  atomic  wave  functions  are 

4 

unambiguously  defined.  Then  it  is  a  standard  procedure  to  eliminate  the 
explicit  occurrence  of  the  radiation  field  from  the  equation  of  motion  (2.4), 
which  yields  the  Einstein  coefficients  for  the  decay  of  |j#m#>.  We  find 
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(2.7) 


where  t  assumes  the  values  -1 ,  0,  1,  corresponding  to  the  three  possible 
helicities  of  a  fluorescent  photon.  The  parameter  functions  b^(d)  embody  the 
surface  modifications  to  the  relaxation  constants,  and  they  are  explicitly 


b|(d)  -  b±J(d)  -  1  "  f  J  du  cos(2(wdu/c  +  9#(ta»,u))) 

-  t  f  du  u2cos(2(wdu/c  +  $  (u,u)))  ,  (2.8) 

*  J  0  p 

bA(d)  ■  bp(d)  -1+|  J  du  (l-u2)cos(2(u>du/c  +  dp(u,u)))  .  (2.9) 


Expressions  (2.8)  and  (2.9)  are  identical  to  those  obtained  by  Philpott,  in  a 
different  way. 

In  the  limit  d  -*■  we  find  by  -  bi  -  1  and  A  -  A  ,  which  implies  that  the 

I  *e  ° 

atom  behaves  as  an  atom  in  free  space  (no  surface),  and  that  the  lifetime  of 
the  state  becomes  independent  of  e#.  Conversely,  due  to  the  presence  of 

the  surface,  every  excited  state  has  a  different  relaxation  constant  A^  .  The 

relation  A^  -  AQ  for  an  atom  in  empty  space  reflects  the  isotropy  of  spontaneous 

e  w 

emission,  or,  the  spherical  symmetry  of  the  system.  In  the  vicinity  of  a 
boundary,  this  symmetry  is  destroyed,  which  gives  every  state  |j#m#>  *  different 
lifetime.  From  Eq.  (2.7)  we  see  that 


A_  -  A 


(2.10) 


which  displays  the  remaining  cylindrical  symmetry  for  rotaions  about  the  z-axis.^ 
By  summing  Eq.  (2.7)  over  the  magnetic  quantum  numbers  m  ,  we  obtain  the  sum 

rule 

JT*i  l  \  -K  +  5*|  •  «•»> 

*  * 

with  Aj_  *  AQbj_  and  A|  ■  AQb|.  Here,  the  left-hand  side  is  the  average  value  of 

A^  over  the  different  substates,  because  there  are  2j@  +  1  excited  states  | 

for  a  given  j#.  If  the  population  of  the  excited  state  is  randomly  distributed 

over  the  various  substates,  with  no  coherence  between  the  states,  then  the  atom 
• 

will  effectively  decay  with  the  average  of  A  .  In  this  respect  the  result 

n 

e 

(2.11)  is  a  quantum  mechanical  justification  for  a  randomization  of  the  dipole 
direction,  which  is  commonly  applied  in  semiclassical  theories  of  relaxation  near 
an  interface.  For  excitation  with  a  polarized  light  source  this  procedure  is  not 
correct,  as  we  shall  show  in  due  course. 

Another  interesting  feature  is  the  dependence  of  the  lifetimes  on  the 
distance  d  between  atom  and  surface.  From  Eqs.  (2.8)  and  (2.9)  we  see  that  b|(d) 
and  bj(d)  are  determined  by  a  superposition  of  waves  with  a  different  value  of  u, 
where  u  is  the  cosine  of  the  angle  of  incidence.  Properties  of  the  reflecting 
medium  are  accounted  for  by  the  phase  shifts  #o(w,u).  From  this  interpretation 
it  follows  that  the  d  dependence  of  A^  is  brought  about  by  simple  interference 
of  plane  waves,  which  is  a  pure  classical  effect.  The  relevance  of  the  phase 
shift  is  most  easily  demonstrated  by  considering  the  limit  d  -*■  0.  In  the  case  of 
a  perfect  conductor  we  then  find 

,  bA(0)  -  2  , 


b|{0)  -  0 


(2.12) 


whereas  the  limit  e  ■*  0  yields 


b||(0)  -  1  ,  bA(0)  -  0  . 


(2.13) 


For  a  perfectly-conducting  solid  the  inverse  lifetime  of  a  parallel  dipole  is 
zero,  but  if  the  conductivity  becomes  smaller,  it  is  the  inverse  lifetime  of  a 
perpendicular  dipole  which  vanishes.  This  reversal  of  behavior  of  parallel  and 
perpendicular  dipoles  as  a  function  of  the  conductivity  illustrates  the  dramatic 
effect  that  a  phase  change  upon  reflection  of  a  plane  wave  can  have  on  the 
optical  properties  of  an  atom,  which  experiences  these  plane-wave  fields.  Figure 
1  gives  more  details  of  the  c  dependence  of  bu(0)  and  bj^(O). 


III.  FLUORESCENCE 

Dynamical  behavior  of  an  atom  in  a  radiation  field  can  be  studied 
experimentally  by  observation  of  the  fluorescence.  Information  on  various 
details,  like  lifetimes,  is  reflected  in  the  different  statistical  properties  of 
the  emitted  radiation.  Therefore,  we  consider  the  (quantum)  radiation  field  in 
more  detail.  In  Eq.  (2.3)  we  gave  the  general  expression  for  the  electric  field 
E(r,t),  and  we  notice  that  its  time  dependence  only  enters  through  the  time 
dependence  of  the  annihilation  operator  a^o(t) .  By  definition,  this  Heisenberg 
operator  obeys  the  time-evolution  equation 


*4  ako(t)  -  [ako(t)  •  H.  +  Hr  +  HiJ 


(3.1) 


If  we  work  out  the  commitator,  differentiate  the  equation  with  respect  to  time 
and  substitute  the  result  in  the  second  time  derivative  of  Eq.  (2.3),  then  it 


•  :'V  • 


■jCh'vv 


8 


follows  that  E(r,t)  obeys  the  wave  equation,  with  the  dipole  moment  as  source 
3 

term.  We  can  solve  this  equation  and  make  an  asymptotic  expansion  for  r  »  |r|  -*• 
•,  which  finally  yields  the  expression  for  the  electric  field  in  the  radiation 
zone 


E(r,t)  -  - — r  V  e  (e  •{ii(t-r/c)e‘i<,,du^c 

4»e  rc2  L  * l 0  ~X~°  " 

o  o 

i( wdu/c+24  (w,u) ) 

+  p*(t-r/c)e  })  +  h.c.  +  free  field  .  (3.2) 

Here,  the  polarization  vectors  correspond  to  the  plane  waves  which  travel  in 
the  direction  of  observation  r,  and  u  is  now  the  cosine  of  the  angle  of 
observation  (angle  between  r  and  the  z-axis).  The  mirror  dipole  moment  p'  is 
defined  as  p'  ■  Pj_  -  Py  in  terms  of  the  separation  p  ■  p^  ♦  py  for  the  atomic 
dipole.  Result  (3.2)  encompasses  all  spectral,  temporal  and  polarization 
properties  of  the  emitted  photons,  including  their  angular  distribution.  We 
notice  that  the  surface-conductivity  effect  is  incorporated  by  the  phase  shift 
of  the  component  with  source  p',  which  originates  from  the  mirror  dipole  in  r  ■  - 
de  below  the  surface. 

IV.  PHOTON  CORRELATIONS 

9 

Recently  we  advocated  that  in  designing  an  experiment  for  the  observation 
of  surface-modified  optical  properties  of  atoms,  one  can  greatly  benefit  from  the 
polarization  properties  of  the  fluorescence,  which  are  induced  by  the  presence  of 
the  substrate.  Especially  the  combination  with  state-selective  excitation  by  an 
intense  polarized  laser  beam,  propagating  in  a  well-defined  direction,  allows  the 
construction  of  a  geometrical  configuration  for  which  only  a  specific  combination 


of  and  determines  the  statistical  properties  of  the  fluorescence.  In  order 

to  illstrate  the  idea,  we  consider  a  j  *  1  •*  J  ■  0  transition.  If  we  take  the 

e  g 

incident  cw  laser  field  to  be  linearly  polarized  and  travelling  in  a  direction 
parallel  to  the  surface,  then  it  follows  from  the  dipole  selection  rules  that 
only  the  state  |l,0>  will  be  populated.  Alternatively,  we  can  assume  a 
circularly-polarized  beam  which  is  incident  normal  to  the  surface.  In  this  case 
only  one  of  the  substates  |l,±l>  will  be  occupied.  Since  |l,0>  and  |l,±l>  decay 
with  an  inverse  lifetime  equal  to  Aj_  and  A||  respectively,  we  can  simply  fix  the 
dipole  direction  with  our  choice  of  external  field. 

Driving  an  atomic  transition  with  a  cw  laser  results  in  a  time- independent 
fluorescent  intensity  I,  and  consequently  the  information  contained  in  this 
quantity  is  minimal.  Therefore,  we  consider  the  stationary  two-photon 
correlation  function  f(t),  which  has  the  significance  of  the  probability  for  the 
detection  of  a  photon  at  time  t  >  0,  after  the  detection  of  a  photon  at  time 
zero,  and  irrespective  of  possible  detections  at  other  times  (for  instance  in 
I0,t]).  In  general  the  function  f(t)  acquires  the  form 

f(t)  «  <E(_)(0)E(’)(t)E(+)(t)E(+)(0)>  ,  (4.1) 

where  E^(t)  is  the  projection  of  E^(r,t)  on  the  polarization  direction  of  a 
polarizer.  Evaluation  of  this  expression  is  quite  involved, ^  so  here  we  shall 
merely  quote  the  result.  We  find 

f(t)  -  I{l-e'3At/4(cos(ftt)+3Asin(ftt)/4Q))  (4.2) 

for  a  strong  laser  exactly  on  resonance  with  the  transition  frequency  u.  Here,  0 

2 

denotes  the  Rabi  frequency  (ft  «  laser  power).  Surface- induced  effects  are 


entirely  incorporated  in  the  Einstein  coefficient  A,  which  equals  Aj_  for  the 
linearly-polarized  case  (An  -  0),  and  Ay  for  a  circularly-polarized  laser  (Am  » 
zl).  The  advantage  of  a  laser-driven  atom  and  observation  of  the  photon 
correlations  is  that  the  desired  relaxation  constant  can  simply  be  selected  by  an 
appropriate  choice  of  external  field. 


V.  MOLECULAR  DECAY  AT  ROUGH  SURFACES 

In  the  previous  sections  perfect  flatness  of  the  surface  has  been  assumed. 
However,  there  are  many  experiments  with  large  admoleucles  (e.g.,  pyrazine)  in 
which  roughness  of  the  surface  unavoidably  exists. Hence,  in  this  section  we 
shall  summarize  some  of  our  recent  work  on  molecular  decay  rates  at  rough 
metallic  surfaces,  where  surface  plasmon  excitation  has  been  found  to  play  a 
dominant  role. 

To  simplify  the  problem,  we  shall  follow  a  classical  phenomenological 

approach  based  on  Maxwell's  theory,  which  has  been  found  to  be  adequate  for 

molecule-surface  distances  (d)  greater  than  a  few  Angstroms.  Then  the  "quantum 

spreads"  of  both  the  surface  and  the  molecule  can  be  neglected.**  Within  this 

approach,  the  molecule  is  modeled  as  a  point  dipole  (p),  satisfying  the  damped- 

oscillator  equation  of  motion,  driven  by  the  surface  field  (E).  The  induced 

12-14 

decay  rate  obtained  in  this  model  can  be  expressed  as 


v  -  v» + !  ^5 


(5.1) 


*R  *  V1 


+  f  Im(G F-K5R)  J 

kJ 


(5.2) 


where  q  is  the  quantum  yield  of  the  emitting  state,  k  *  u/c  equals  the  emission 
wave  number,  and  the  sub( super )scripts  R,  F  and  o  stand  for  cases  with  a  rough 


surface,  flat  surface  and  free  molecule  respectively.  The  G-functions  in  Eqs. 
(5.1)  and  (5.2)  are  defined  as  the  surface  field  acting  on  the  molecule  per  unit 
dipole  moment  and  as  a  function  of  the  emission  frequency 


G(w)  -  ^  .  (5.3) 

Note  that  in  writing  Eq.  (5.2)  we  have  followed  an  approach  which  distinguishes 

between  the  contributions  to  the  surface- induced  decay  rates  from  the  flat 

13  14 

boundary  and  from  the  roughness,  respectively.  ' 

The  remaining  task  in  this  approach  is  to  calculate  the  G(u)-functions.  In 

general,  there  are  two  ways  of  doing  this,  namely,  the  image  method  and  the 

energy-transfer  theory.  In  the  image  theory  (IT),  one  regards  the  surface  field 

E(u)  as  given  by  the  image  field,  induced  by  the  molecular  dipole,  which  is 

obtained  from  solving  the  electrostatic  Poisson  equation.  For  example,  in  the 

case  of  a  flat  surface  and  y  located  at  (0,0,d),  oriented  perpendicular  to  the 

14 

surface,  one  obtains 


(<*>) 


e-1 

e+1 


(5.4) 


where  3  *  kd  is  the  reduced  distance  and  c(w)  ■  e^(w)  +  ie2(w)  is  the  complex 
bulk  dielectric  constant  of  the  substrate  material. 

The  other  approach,  energy-transfer  theory  (ET),  is  a  more  exact 
electrodynamical  treatment  in  which  one  has  to  solve  the  full  set  of  Maxwell's 
equations  (or  the  Helmholtz  wave  eqution)  by  regarding  the  dipole  emission  as  the 
source  (incident)  field,  and  by  matching  boundary  conditions  at  the  surface. 
Then  E(w)  will  be  given  by  the  reflected  field  from  the  surface,  evaluated  at  the 
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molecular  sit*.  For  th«  can  of  s  fist  surface,  this  problem  has  been  solved  by 

12 

Chance,  Frock  and  Silbey  (CPS)  by  applying  the  Soonerfeld  theory  for  a 
radiating  dipole  antenna  above  the  earth's  surface.  Thus  for  the  above  example 
the  CPS  theory  gives 


cJjU)  -  -k3  fdu  RpU.ilj) 


-21,3  3 
1  u 


(5.5) 


where  ■  -i(l-u3)^. 

It  is  obvious  from  the  above  examples  that  IT  is  in  general  an  approximate 

but  much  simpler  approach  than  ET.  It  has  been  a  wide  belief13  that  in  the 

short -distance  or  long-wavelength  limit  (3«1),  which  holds  for  most 

2  3  16 

experimental  conditions  (d  <  10  A  and  X  -  10  A),  IT  should  be  accurate  enough 

to  describe  the  phenomenon.  Moreover,  although  both  IT  and  ET  are  available  for 

the  case  of  flat  surfaces,  only  IT  is  available  so  far  in  the  literature  when 
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surface  roughness  cannot  be  neglected.  * 

In  a  recent  analysis  of  the  limit  of  IT,  however,  it  was  pointed  out  that 

the  condition  d  «  A  is  not  sufficient  to  ensure  IT  to  be  valid. 1^  The 

additional  condition  found  for  IT  to  hold  is  d  «  6,  with  6  being  the  skin  depth 
of  the  substrate.  In  particular,  for  a  highly-conducting  substrate  such  as 
silver  and  for  the  above  example,  IT  has  been  found  to  break  down  appreciably  for 
d  >  0.01A.  Roughly  speaking,  the  physical  origin  for  such  a  failure  of  IT  stems 
from  the  fact  that  the  Helmholtz  wave  equation  in  a  conducting  medium  does  not 
reduce  to  the  static  Poisson  equation  in  the  limit  X  -*■  •,  if  the  conductivity  of 
the  medium  becomes  infinite. 

It  is  therefore  necessary  to  develop  a  dynamical  theory  (ET)  for  a  complete 

description  of  the  decay  rates  of  admolecules  at  rough  surfaces.  Recently,  an 

18  12 
attempt  in  this  direction  was  made  by  combining  the  CPS  theory  and  the 


integral -equation  formalism  of  Maxwell's  equations  for  rough  boundaries, 
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established  chiefly  by  Maradudin  and  Mills  (MM)  and  Agarval.  Let  us  describe 
the  roughness  by  a  profile  function  C(x,y).  The  main  idea  of  the  MM  theory  is  to 
regard  £  as  a  source  for  the  homogeneous  Hemholtz  equation  in  the  case  of  a 
perfectly-f lat  boundary.  One  is  then  able  to  obtain  perturbative  solutions  in 
orders  of  C  for  the  reflected  field  E(w),  which  is  the  relevant  surface  field  for 
the  calculation  of  the  G(u)-functions.  Then  E(w)  is  expressible  as  a  combination 
of  the  contributions  from  the  flat  boundary  and  the  roughness  as 

E  -  E*  +  E*  ,  (5.6) 


and  to  order  C  one  obtains  for  the  p-th  component  of 
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Ej(rju.)  -  -  -^3  [e(«)-l]  f  d2k|  e1’*  £<k(-kjj0)) 

*  j  d*'  d^kyulz*'  )«(*')  E^^(k|^w|s* )  .  (5.7) 

Here,  E^  is  the  total  field  for  a  perfectly  flat  boundary,  and  £  and  d^v  are 

the  Fourier  transforms  of  the  profile  function  and  the  Green  tensors, 

respectively,  as  given  in  Ref.  19.  The  evaluation  of  the  6  functions  for 

discontinuous  fields  across  a  boundary  is  understood  to  follow  MM  together  with 
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Agarwal's  prescription.  ' 

To  illustrate  this  dynamical  theory,  we  consider  the  simple  case  of  a 
perpendicular  dipole  located  at  (0,0, d)  above  a  shallow  sinusoidal  grating. 
After  application  of  the  CPS  theory  for  the  calculation  of  E^,  we  obtain  to 
lowest  order  in  C 


,  (5. 


2  ik2d 
kue 


^<d:“)  ■  ■iuk3  ^  1  dN  «*!>  C du  “‘V  i  •  l‘a 


where 


kx  -  -[e(w)k2-kjjj*  , 


k2  " 


(k2-kjj)*  ,  k2  >  kjj 
i(kjj-k2)*  ,k2  <  kjj 


(5.9) 


For  C(x,y)  ■  Q  e*^*  we  obtain 
o 


£<k||)  -  (2s)2  Cc  «(Q-kj|)  , 


(5.10) 


where 


9-Q?x  • 


(5.11) 


Equation  (5.8)  finally  yield*18 


2  ik2d 


<4  -  i  ■  -i*3  f0  J“  <“V'  *“ia  • 


(5.12) 


where  kj  and  k2  are  now  defined  as  in  Eq.  (5.9)  with  ky  replaced  by  Q. 


IS 


Previously,  we  have  also  solved  the  sane  problem  by  the  image  method  through 

the  application  of  the  image  theory  for  a  rough  boundary,  established  by  Rahman 
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and  Maradudin.  The  result  thus  obtained  can  be  expressed  as 

G?t  ■  ~ C  l  fdu  f°dv  (efg+h)exp[-(f+g)dj  ,  (5.13) 

IT  *  (e+1)2  J0  J0 

where  f,  g  and  h  are  functions  of  u  and  v,  given  by 
f(u,v)  -  [(u+S)2  +  v2]*  , 

g(u, v)  -  [(u-^)2  +  v2]^  ,  (5.14) 

h(u,v)  ■  u2  +  v2  -  . 

Hence  for  this  problem,  according  to  Eq.  (5.6),  Eqs.  (5.4)  and  (5.13) 
provide  a  complete  description  of  A  within  IT,  whereas  Eqs.  (5.5)  and  (5.12) 
provide  a  complete  dynamical  theory  (ET)  for  A.  We  have  carried  out  detailed 
numerical  studies  of  these  results  for  different  metal  substrates  such  as  Ag,  Au, 
Cu  and  Ni,  where  both  w  and  d  were  varied.  For  illustrative  purposes,  we  have 
reproduced  some  results  in  Fig.  2.  In  general,  we  arrive  at  the  following 
conclusions: IT  can  be  very  inaccurate  for  both  flat  and  rough 
surfaces  in  the  case  of  a  highly-conducting  substrate  or  for  large  molecule- 
surface  distances;  (ii)  the  surface-plasmon  effect  can  cause  serious  damping  when 
c.  is  close  to  the  resonance  frequency;  (iii)  surface  roughness  can  either  enhance 
cr  diminish  the  flat-surface  value  for  this  induced  decay,  depending  on  u,  d  and 


the  orientation  of  the  admolecule. 


Although  we  have  most  of  the  time  considered  the  sinusoidal  grating  as  a 

prototype  of  rough  surfaces,  other  kinds  of  roughness  have  also  been  considered 

in  the  literature.  As  an  example.  Arias,  Aravind  and  Metiu^  have  constructed  an 

IT  for  admolecules  at  surfaces  with  randomly-distributed  (Gaussian)  roughness. 

To  extend  ET  to  this  case,  one  has  to  include  contributions  to  E(u)  which  are 

18 

second  order  in  C»  since  the  first-order  contribution  will  vanish  in  this  case. 
Nevertheless,  we  expect  that  similar  comparisons  will  be  obtained  as  in  Fig.  2. 

VI.  CONCLUSIONS 

In  this  paper,  we  have  summarized  some  of  our  recent  results  on  the 

spectroscopic  properties  for  atoms  as  well  as  molecules  in  the  vicinity  of  a 

dielectric  surface.  For  the  case  of  adatoms  near  a  flat  surface,  a  fully  quantum 

electrodynamical  (QED)  treatment  has  been  followed  where  the  atomic  lifetimes, 

polarization  of  the  atomic  dipole  and  photon  correlation  phenomena  have  been 

analyzed.  For  the  case  of  a  molecular  dipole  in  the  vicinity  of  a  rough  metallic 

surface,  a  simpler  classical  phenomenological  approach  has  been  followed  to 

establish  a  dynamical  theory  of  the  surface-induced  decay  rates.  It  is  found 

that  the  hitherto  available  static  (image)  theory  has  very  limited  validity. 

Future  directions  in  this  research  may  involve  a  more  exact  QED  treatment  for  the 

case  of  rough  surfaces  (along  the  same  lines  as  we  did  for  flat  surfaces),  as 

well  as  possible  nonlinear  surface  spectroscopy  due  to  a  high- intensity  laser 

field  or  special  surface  properties.  To  this  later  aspect,  we  have  recently 

looked  into  the  possibility  of  applying  the  technique  of  optical  phase 

22 

conjugation  by  a  surface  to  study  possible  novel  spectroscopic  phenomena  for  a 
system  of  molecules  in  the  vicinity  of  such  surfaces. 
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FIGURE  CAPTIONS 


Fit.  1.  Plot  of  the  parallel  (curve  a)  and  perpendicular  (curve  b)  inverse 

lifetimes  as  a  function  of  the  relative  permittivity  e,  and  in  the  limit  of  a 

zero  atom-surface  distance.  For  c  ■*  -•  the  curves  approach  the  values  b||(0)  ■  0 
|  and  bA(0)  ■  2,  corresponding  to  a  perfect  conductor.  We  notice  that  b|(0)  is 

i  almost  independent  of  c  and  approximately  equal  to  zero  over  the  entire  c  range. 

! 

Only  for  e-  0  does  it  attain  a  finite  value.  From  curve  b  we  see  that  bj^O) 

gradually  changes  as  a  function  of  c(  and  that  it  never  reaches  its  perfect- 

) 

i  conductor  limit  for  any  reasonable  value  of  c. 

Fit.  2.  Comparison  between  the  energy  transfer  theory  (ET,  solid  curves)  and  the 
image  theory  (IT,  dotted  curves)  for  decay  rates  of  admolecules.  Curves  a  and  b 
are  the  values  of  IaG  for  an  Ag  grating  substrate  ((q  ■  0.7A,  Q  ■  0.01  A  )  at  u 

S  F 

i  >2.5  eV,  and  curves  c  and  d  are  the  values  of  ImG  for  a  flat  Ni  surface  at  u  ■ 

3.3  eV.  The  unit  of  G  is  A  .  The  role  of  the  substrate  conductivity  is  evident 
from  these  graphs. 
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